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Abstract 

The bulk and shear viscosities of a pion gas is obtained by solving the relativistic transport 
equation in the Chapman-Enskog approximation. In-medium effects are introduced in the tttt 
cross-section through one- loop self-energies in the propagator of the exchanged p and a mesons. 
The effect of early chemical freeze-out in heavy ion collisions is implemented through a temperature 
dependent pion chemical potential. These are found to affect the temperature dependence of the 
bulk and shear viscosities significantly. 
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I. INTRODUCTION 



The study of transport coefficients and non-equilibrium dynamics in general has assumed 
a great significance in recent times. The impetus in this direction has been provided by 
the recent results from the Relativistic Heavy Ion Collider (RHIC), in particular the elliptic 
fiow which seems to be rather well described by nearly ideal hydrodynamics with a very low 
value of Tj/s, close to the quantum bound l/An [ij], rj and s being the coefficient of shear 
viscosity and entropy density res|3ectively. This led to the description of quark gluon plasma 
as the most perfect fiuid known [2]. This bound is however found to be violated, e.g. in the 
case of a strongly coupled anisotropic plasma (sl. 

Of the two viscosities used to parametrize the leading order corrections to the stress tensor 
the shear viscosity rj is more commonly discussed, the magnitude of C, the bulk viscosity 
usually being much smaller in comparison. The latter actually vanishes as in the case of 
systems with conformal invariance e.g. a gas of non-interacting highly relativistic particles. 
Breaking of conformal invariance of QCD due to quantum effects shows up through non-zero 
values of C, which is related to the trace of the energy momentum tensor. An estimate of this 
can be obtained from the interaction measure determined from lattice calculations ^ . The 
behaviour of C and 77 as a function of temperature is particularly relevant in the context of 
non-ideal hydrodynamic simulations of heavy ion collisions. Whereas rj/s as a function of 
T is expected to go through a minimum 2!] at or near the critical value for crossover from 
the partonic phase, it is believed that C/s may be large |^ or even diverging close to this 
value. These have been studied in the linear sigma model in the large- iV limit jo, Q] and for 
a massless pion gas in js]. 

The viscosities of a pion gas have received some attention in recent times. In the dia- 
grammatic approach, one uses the Kubo formula which relates the transport coefficients to 
retarded two-point functions j^, 1^. The shear viscosity of a pion gas in this approach was 



evaluated in 11 



12l |. However, the kinetic theory approach being computationally more 



efficient 13j], has been mostly used to obtain the viscous coefficients. The vrvr cross-section 
is a crucial dynamical input in these calculations. Scattering amplitudes evaluated using 
chiral perturbation theory fl^. \\\ to lowest order have been used in 16, 17| and unitariza- 



tion improved estimates of the amplitudes were used in [18| to evaluate the shear viscosity. 
Phenomenological scattering cross-section using experimental phase shifts have been used 
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in 171, Il9l42l|. While in 22| the effect of number changing processes on the bulk viscosity 



of a pion gas has been studied, in 23| unitarized chiral perturbation theory was used to 
demonstrate the breaking of conformal symmetry by the pion mass. 

Our aim in this work is to construct the tttt cross-section in the medium and study its effect 
on the temperature dependence of the viscous coefficients. For this purpose we employ an 
effective Lagrangian approach which incorporates p and a meson exchange in vrvr scattering. 
A motivating factor is the role of the p pole in vrvr scattering in preserving the quantum 



bound on T]/s for a pion gas as demonstrated through use of the phenomenological 20|, 124 1 



and unitarized 18| cross-section in ^5j. Medium effects are then introduced in the p and a 
propagators through one-loop self-energy diagrams. Now, the hadronic matter produced in 
highly relativistic heavy ion collisions is known to undergo early chemical freeze-out. Number 
changing (inelastic) processes having much larger relaxation times go out of equilibrium at 
this point and a temperature dependent chemical potential results for each species so as 
to conserve the number corresponding to the measured particle ratios. We thus evaluate 
the bulk and shear viscosity of a pion gas below the crossover temperature in heavy ion 
collisions considering only elastic vrvr scattering in the medium with a temperature dependent 
pion chemical potential. In the process we extend our estimation of the shear viscosity at 



vanishing chemical potential 26| where a significant effect of the medium on the temperature 
dependence was observed. 

In the next section, we briefly recapitulate the expressions for the viscosities obtained by 
solving the transport equation in the Chapman-Enskog approximation. In section III we 
describe the vrvr cross-section in the medium. Numerical results are provided in section IV 
followed by summary and conclusions in section V. Various details concerning the solution 
of the transport equation are provided in Appendices A, B and C. 

II. THE BULK AND SHEAR VISCOSITY IN THE CHAPMAN-ENSKOG AP- 
PROXIMATION 

The evolution of the phase space distribution of the pions is governed by the equation 

p''d,f{x,p) = C[f] (1) 
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where C[f] is the colhsion integral. For binary elastic collisions p + k p' + k' which we 



consider, this is given by 



2l| 



C[f] = JdTkdTp, dTk'[fix,p')f{x,k'){l + fix,p)}{l + f{x,k)} 

-fix,p)f{x,k){l + fix,p')}{l + f{x,k')}] W (2) 



where the interaction rate, 

s da 



and dVq = i^^yi^ with E = JfP + m^. The 1/2 factor comes from the indistinguishability 



(27r)3£; 

of the initial state pions. 

For small deviation from local equilibrium we write, in the first Chapman- Enskog ap- 
proximation 

f{x,p) = f'\x,p) + Sfix,p), 5f{x,p) = f^'\x,p)[l + f^'\x,pMx,p) (3) 
where the equilibrium distribution function is given by 



f(^\x,p) 



(4) 



with T(x), u^{x) and ^■,^{x) representing the local temperature, flow velocity and pion chem- 
ical potential respectively. Note that the metric diag{l, —1, —1, —1) is used. Also, we take 
Ufj,u^ = 1 where = (1, 0) in the local rest frame. Putting ([3]) in ([T]) the deviation function 
(p{x,p) is seen to satisfy 

p^d,f'\x,p) = -£[0] (5) 
where the linearized collision term 



m = f'\x,p) J dr, dr,> dr,,f'\x,k){i + f'\x,p')}{i + f'\x,k')} 

[(f){x,p) + k) — (j){x,p') — k')] W . (6) 

The form of as given in (jlj) is used on the left side of ([5]) and the time derivatives 

are eliminated with the help of equilibrium thermodynamic laws. As detailed in Appendix-B 
this leads us to the equation 

[Qd,u'' + p,A>^''{p ■ u - h){T-^d,T - Du,) - {p,Pu){d^u'')]f^'^\l + /(°)) = -T£[0] (7) 
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where D = uf'd^, = A^^d", A^,^ = g^^ - m^m^ and the notation {tf"") = [^{A^"^A''^ + 
A'^^A^'^) — lA^'^A^^Jto,^ indicates a space-hke symmetric and traceless form of the tensor 
f^". In this equation 

Q = -Iml + (p . uf{^ - 7'} + P ■ «{(7" - l)h - Y'T} (8) 



where 



{SI/ Sly - {Sl/S',)^ + 3z-'S',Sl/{Siy + 2z-^Sl/Sl - 

III S\IS\-\-hz ^»S'3/»S'2 S'^Sk^ I {S}^'^ /I 1 \ 

^ ^ (^27^2^)' - {SllSlf + 3^-150537(51)2 + 2z-'SllS\ - ^ ^ 

with 2; = m^r/T and /i = mTr^l/Sj. The terms 5" are integrals over Bose functions and 

are defined in Appendix-A. The left hand side of ([5]) is thus expressed in terms of the 

thermodynamic forces d ■ m, A^^ {T^^dyT — Duy) and {d^u") which have different tensorial 

ranks, representing a scalar, vector and tensor respectively. In order to be a solution of this 

equation (p must also be a linear combination of the corresponding thermodynamic forces. 

It is typical to take as 

(f) = Ad-u + B^A^^iT-^d.T - Du,) - C^^id'^u") (12) 

which on substitution into ([7]) and comparing coefficients of the (independent) thermody- 
namic forces on both sides, yields the set of equations 

£[A] = -Q/W(p){l + /(°)(p)}/T (13) 

Cn.] = -{p,Pu)f^'\p){l + f^'\p)}/T (14) 

ignoring the equation for 5^ which is related to thermal conductivity. These integral equa- 
tions are to be solved to get the coefficients A and C^y. It now remains to link these to 
the viscous coefficients ( and t]. This is achieved by means of the dissipative part of the 
energy-momentum tensor resulting from the use of the non-equilibrium distribution function 
m in 

T^'' = J dTp p^p^fip) = T'^^^") + AT'^^ (15) 

where 

AT^" = J rfrp/(°)(l + /^°))C„;3(pV)(5''m") + / dTpf^'^\l + /(°))gAA'^"9 ■ u (16) 



Again, for a small deviation (f){x,p), close to equilibrium, so that only firs , 
contribute, the dissipative tensor can be generally expressed in the form 



order derivatives 
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(17) 



Comparing, we obtain the expressions of shear and bulk viscosity, 
and 



(18) 



C = -JdT,QAf^'\p){l + f'\p)}. (19) 

The coefficients A and C^jy are perturbatively obtained from fll3l) and (1141) by expanding in 
terms of orthogonal polynomials which reduces the integral equations to algebraic ones. This 

281]. A brief account is 



elaborate procedure using Laguerre polynomials, is described in |21 
provided in Appendix C. Finally, the first approximation to shear and bulk viscosity come 
out as 
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and 



where 



c = 

022 



7o 
Coo 



-10-"3(^) 



lQ{I,{z) + h{z) + -h{z)} 



and 



Q!2 



^22 



qO o cl cl 

- ^-') + + -|t){(i - 7")| + rz-'} 



2 ^3^52 
4 



5° 



5*2 Z 5*2 

I On 



2z'h(z) . 



(20) 
(21) 



(22) 



(23) 



The integrals Ia{z) appearing in the above expressions are defined as 

.4 I 



Uz) 



dx sinh^" x I dO sin 9 — g — — — -r^-^^ -r 

Jo [e^ - l){e^ - l){e^ - l){e" - 1) 



/ dip cosh^ '?/' sinh''' / dQ sm&- — (ilj,Q) / d(j) 
Jo Jo 2d\l Jo 



^2z cosh 1/1 cosh x 



M^{9,Q) (24) 
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FIG. 1: The tttt cross-section as a function of centre of mass energy. The dotted hne indicates 
the cross-section obtained using eq. ()29p which agrees well with the experimental values (eq. (j26p ) 
shown by filled triangles. The dashed and solid lines depict the in- medium cross-section for vrvr and 
irh (= 7r,cj, /ii, oi) loops respectively in the p and a self-energies evaluated at r=160 MeV with 
temperature dependent pion chemical potential. 

in which the functions Ma{0, 0) represent 

Mi{e,e) = 1 -cos^e , 

Ms (0, 6) = cos^ 9 + cos^ 9' -2 cos 9 cos 9' cos 9 , 

M3{9,Q) = [cos^ ^ - cos^ . (25) 

Note that the differential cross-section which appears in the denominator is the dynamical 
input in the expressions for t] and (. It is this quantity we turn to in the next section. 



III. THE vrvr CROSS-SECTION WITH MEDIUM EFFECTS 



The strong interaction dynamics of the pions enters the collision integrals through the 
cross-section. In Fig. [T] we show the vrvr cross-section as a function of the centre of mass 
energy of scattering. The different curves are explained below. The filled triangles referred 



to as experiment is a widely used resonance saturation parametrization 
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FIG. 2: (Color online) The tttt phaseshifts (I27p as a function of E compared to Roy equation 



solutions given in Ref. 



30| 



and isovector phase shifts obtained from various empirical data involving the tttt system. 
The isospin averaged differential cross-section is given by 



dcr{s) 



4 



^ sin^ ^0 + - sin^ 6^ + - - 9 sin^ 6l cos^ 9 



9 



9 



(26) 



where 



TT 



+ arctan 



TT 

Sl = — + arctan 



' E — rricr 
E — nip" 



The widths are given by = 2.06J9 and Tp = OMbp 



with vHn 



(27) 



5.8m^ and 



VTLr 



5.53m^ . As seen in Fig. [2] these phase shifts agree quite well with those obtained 



from solutions of the Roy equations as given in [30| . The bands bordered by the dotted lines 
represent the uncertainties in the solution. The experimentally measured phase shifts (not 
shown) have error bars [30] which are not reflected in the parametrizations (1271) plotted in 

Fig.m 

Our objective now is to set up a dynamical model which agrees reasonably with the above 
parametrization in vacuum and is amenable to the incorporation of medium effects using 
many body techniques. We thus evaluate the tttt cross-section involving p and a meson 
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exchange processes using the interaction Lagrangian 



(28) 



where Qp = 6.05 and = 2.5. In the matrix elements corresponding to s-channel p and a 
exchange diagrams which appear for total isospin 1=1 and respectively, we introduce a 
decay width in the corresponding propagator. We get 26| . 



Mi=o 
Mi=i 

Ml=2 



2 2 



2 2 



2 2 



1 



1 



t — rri^ 
1 



u — mi 



t — mi 

,2 



t 



+ 



mt 



u — mt 



9o 



+ 9; 



u — m^ 
2{t - u) 



^9t 



s — u s 

+ 



s — + imoT 



+ 



t — m^p 
t - s 



u 



+ 



mi 



t-s 
u — mi 



u — mi 



u — mi 



u — s 
t — 



(29) 



The differential cross-section is then obtained from ^ 



|A^p/647r s where the isospin 



averaged amplitude is given by = 1 1] (2/ + 1)|A^/P. 

Considering the isoscalar and isovector contributions which involves s-channel a and p me- 
son exchange respectively, the integrated cross-section, shown by the dotted line (indicated 
by 'vacuum') in Fig. [T]is seen to agree reasonably well with the experimental cross-section 
up to a centre of mass energy of about 1 GeV beyond which the theoretical estimate gives 
higher values. We hence use the experimental cross-section beyond this energy. Inclusion 
of the non-resonant 1 = 2 contribution results in an overestimation and is ignored in this 
normalization of the model to experimental data. It is essential at this point to emphasize 
that the present approach of introducing finite decay widths for the exchanged resonances 
explicitly is not consistent with the well-known chiral low energy theorems concerning the 
amplitude for vrvr scattering 3l|] as well as the scattering lengths. 

We now construct the in-medium cross-section by introducing the effective propagator 
for the p and a mesons in the above expressions for the matrix elements. We first discuss 
the case of the p followed by the relatively simpler case of the a meson. 

The in-medium propagator of the p is obtained in terms of the self-energy by solving the 
Dyson equation and is given by 



flU 



(30) 



where D^^^ is the vacuum propagator for the p meson and 11'^^ is the self energy function 
obtained from one-loop diagrams shown in Fig. [31 The standard procedure 32| to solve this 



FIG. 3: Self-energy diagrams where h stands for 7r,u},hi,ai mesons. 



equation in the medium is to decompose the self-energy into transverse and longitudinal 
components. For the case at hand the difference between these components is found to be 
small and is hence ignored. We work with the polarization averaged self-energy function 
defined as 



where 



n 



n 



2 • - 

The in-medium propagator is then written as 



(31) 



(32) 



(33) 



-m^- Ren(go, q) + 2lmn(go, q) 

The scattering, decay and regeneration processes which cause a gain or loss of p mesons in 
the medium are responsible for the imaginary part of its self-energy. The real part on the 
other hand modifies the position of the pole of the spectral function. 

In the real-time formulation of thermal field theory the self-energy assumes a 2x2 matrix 
structure of which the 11-component is given by 



nil(g) = ^ j 



(2vr) 



(34) 



where D^^ is the 11-component of the scalar propagator given by D^^{k) = A{k) + 
2t: if^^\k) 5 {k"^ — w?). It turns out that the real and imaginary parts of the self-energy 
function which appear in eq. ( l33l) can be obtained in terms of the 11-component through 



the relations (321 . 



33] 



Ren 



Ren 



11 



imn 



e(go)tanh(/3go/2)Imn 



11 



(35) 



Tensor structures associated with the two vertices and the vector propagator are included 

3J] where the interactions were taken from chiral perturbation 



in Nfj^i, and are available in 



10 



theory. It is easy to perform the integral over ko using suitable contours to obtain 
d^k 1 



go - Wtt + Wfe + «^?e(go) 



^ /W(a;.)iVr-/^°)K)iV3"'^ ^ -/(°)K)iVr - (1 + /(°)K))iV; 



qo + uj^ - ujh + imiQo) qo + + Uh + iv^iqo) 

where f^'^^co) = ^(,j-^1)/t_i is the Bose distribution function with arguments = 
and Uh = \l {q — kY + Tn\. Note that this expression is a generalized form for the in-medium 



(36) 



k'^ + mi 



self-energy obtained by Weldon 



35| . The subscript i{= 1, ..4) on N^'^ in fl36l) correspond to 



its values for k^ 



-Ut^, qo — ijJh,qo + (^h respectively. It is easy to read off the real and 



imaginary parts from ( !36l) . The angular integration can be carried out using the (5-functions 
in each of the four terms in the imaginary part which define the kinematically allowed regions 
in go and q where scattering, decay and regeneration processes occur in the medium leading 
to the loss or gain of p mesons 3J]. The vector mesons u, hi and ai which appear in the loop 
have negative G-parity and have substantial Sir and pvr decay widths [36]. The (polarization 
averaged) self-energies containing these unstable particles in the loop graphs have thus been 
folded with their spectral functions, 

1 



dM^-lm 

(mfe-2rft)2 TT 



M^-ml + iMTh{M) 



n(g,M) 



(37) 



with Nh 



(mh+2rh)2 I 

rfM -Im 

(m;,-2rh)2 vr 



The contributions from the loops 



-ml + iMTh{M) 

with heavy mesons (the nh loops) may then be considered as a multi-pion contribution to 
the p self-energy. 

The medium effect on propagation of the a meson is estimated analogously as above. 
The effective propagator in this case is given by 



D{qo,q) 



-1 



q^ - ml - Ren(go, q) + iImn(go, q) 
Following the steps outlined above the expression for the self-energy of the a is given by 



(38) 



n(go,g) 



N 



+ 



d^k 



1 + /W(a;.) + /(o)(a;;) , /W(u;;) - /^(a;.) 



+ 



go - - ir]e{qo) go - + + «77e(go) 



/W(a;.)-/W(a;;) i + /(o)(^^) + /(o)(^;) 



go + - + «^?e(go) go + + + i?7e(go) 



(39) 



where uj'^ = y {q — ky + m^. The imaginary part for the kinematic region of our interest in 
this case receives contribution only from the first term which essentially describes the decay 
of the a into two pions minus the reverse process of formation. 
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FIG. 4: The pion chemical potential as a function of temperature 39l |. 



The in-medium cross-section is now obtained by using the full p and a propagators 
given by (!33|) and (!38l) respectively in place of the vacuum propagators in the scattering 
amplitudes. The long dashed line in Fig. [1] shows a suppression of the peak when only the 
vrvr loop in the p and a self-energies are considered. This effect is magnified when the irh 
loops containing heavier mesons in the p self-energy are taken into account and is depicted 
by the solid line indicated by multi-pion. This is also accompanied by a small shift in the 
peak position. Extension to the case of finite baryon density can be done using the spectral 
function computed in where an extensive list of baryon (and anti-baryon) loops are 
considered along with the mesons. A similar modification of the vrvr cross-section for a hot 



and dense system was seen also in 



We end this section with a discussion of the pion chemical potential. As mentioned in the 
introduction, in heavy ion collisions pions get out of chemical equilibrium early at T ~ 170 
MeV and a corresponding chemical potential starts building up with decrease in temperature. 
The kinetics of the gas is then dominated by elastic collisions including resonance formation 
such as vrvr ^ p etc. At still lower temperature, T ~ 100 MeV elastic collisions become rarer 
and the momentum distribution gets frozen resulting in kinetic freeze-out. This scenario is 
quite compatible with the treatment of medium modification of the vrvr cross-section being 
employed in this work where the vrvr interaction is mediated by p and a exchange and the 
subsequent propagation of these mesons are modified by two-pion and effective multi-pion 



fluctuations. We take the temperature dependent pion chemical potential from Ref. 



39| 
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FIG. 5: (Color online)The bulk viscosity in various scenarios as a function of T. The upper (with 
circles), middle (with triangles) and lower sets of curves correspond to /i^r = 85 MeV, /Utt = fJ-niT) 
and fijr = respectively. In each set the solid line represents use of vacuum cross-section, the 
dotted line for in-medium modification due to pion loop and the dashed line for loops with heavy 
particles in addition. 



which implements the formalism described in 40| and reproduces the slope of the transverse 



momentum spectra of identified hadrons observed in experiments. Here, by fixing the ratio 
s/n where s is the entropy density and n the number density to the value at chemical freeze- 
out where yU^r = 0, one can go down in temperature up to the kinetic freeze-out by increasing 
the pion chemical potential. This provides the temperature dependence leading to fin{T) 
which is shown in Fig. |H In this partial chemical equilibrium scenario of 40| the chemical 
potentials of the heavy mesons are determined from elementary processes. The u chemical 
potential e.g. is given by yU^^ = 3 x 0.88/2,^, as a consequence of the processes u ^ tttttt 



occurring in the medium. The branching ratios are taken from 36 1. 



IV. RESULTS 



We begin with the results for bulk viscosity C as a function of temperature T. In Fig. [5] 
the three sets of curves correspond to different values of the pion chemical potential. The 
uppermost set of curves (with circles) show the bulk viscosity calculated with a pion chemical 
potential //^ ~ 85 MeV. The corresponding curves in the lowermost set are evaluated with 
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FIG. 6: The entropy density of an interacting pion gas as a function of T for different values of 
the pion chemical potential. 

/iir = 0. These values are representative of the kinetic and chemical freeze-out in heavy 
ion collisions respectively. The solid line in the lowermost set represents the case where the 
vacuum cross-section given by eq. (l29l) is used and agrees with the estimate in j21|. The set of 
curves with triangles depicts the situation when fi^^ is a (decreasing) function of temperature 
as shown in Fig. |H This resembles the situation encountered in the later stages of heavy ion 
collisions and interpolates between the results with the constant values of the pion chemical 
potential discussed before. The three curves in each set show the effect of medium on the 
TTTT cross-section. The short-dashed lines in each of the sets depict medium effects for pion 
loops in the p propagator and the long dashed lines correspond to the situation when the 
heavy mesons are included i.e. for nh loops where h = 7r,uj,hi,ai. The clear separation 
between the curves in each set displays a significant effect brought about by the medium 
dependence of the cross-section. A large dependence on the pion chemical potential is also 
inferred since the three sets of curves appear nicely separated. 

Viscosities for relativistic fluids are generally expressed in terms of a dimensionless ratio 
obtained by dividing with the entropy density. The latter is obtained from the thermody- 
namic relation 

Ts = e + P -rifi^ . (40) 
For a free pion gas, we get on using the relations for the energy density e, pressure P and 
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FIG. 7: as a function of T for different vrvr cross-section. The temperature dependent pion 
chemical potential has been used in all cases. 



number density n from Appendix-A, 



27r2 



(41) 



Here g,, = 3, z = mT,/T and the functions S"{z) are given in Appendix-A . Interactions 
between pions lead to corrections to this formula. To O(T^) this has been calculated for 
finite pion chemical potential in 4l|] using chiral perturbation theory to give 



As 



3m 



(42) 



where /„ 



given m 
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93 MeV. It is easily verified that this expression reduces for /i^ = to those 



42]. This correction is ~ 1 — 2% for values of ^-^ and T considered here. In 
Fig. [6] the entropy density of an interacting pion gas as a function of temperature is shown 
for three values of the pion chemical potential as discussed in the context of Fig. [51 

In Fig. [7] we show C/s as a function of T using the temperature dependent pion chemical 
potential. The medium dependence is clearly observed when we compare the results obtained 
with the vacuum cross-section with the ones where the a and p propagation is modified due 



to vrvr and vr/i (multi-piori) 



observed also in [22] and |43 |. 



oops. The decreasing trend with increasing temperature was 



We now turn to the shear viscosity. Here we extend the work in Ref. 26|] to include the 
effect of the pion chemical potential. Shown in Fig. [8] is the shear viscosity as a function of 
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FIG. 8: rj as a function of T for different values of the tttt cross-section. The temperature dependent 
pion chemical potential has been used in all cases. 
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FIG. 9: (Color online)ry/s as a function of T. The upper set of curves with circles correspond to 
/iir = and the lower set corresponds to fij^ = Hn{T). 

T where the results with tttt and nh loops are contrasted with the case where the vacuum 
cross-section is used. The result with the vacuum cross-section agrees with 



and 



2M for 



yU^r = 0. A noticeable medium effect is observed as indicated by the short and long-dashed 
lines. 

Finally, in Fig. |9]we show results for rj/s. For /i^ = the upper set of curves with filled 
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circles show the usual decreasing trend as seen, for example in 4J]. This trend is reversed 
when fin{T) is used and rj/s increases with T in contrast to (/s which decreases. The values 
in all cases remain well above l/47r. 



V. SUMMARY AND OUTLOOK 

In this work we have evaluated the viscosities of a pion gas by solving the Uehling- 
Uhlenbeck transport equation in the Chapman-Enskog approximation to leading order with 
an aim to study the effect of a medium dependent cross-section. The vrvr cross-section 
which goes as an input to these calculations is evaluated from a and p exchange processes. 
Spectral modifications of the p and cr propagators implemented through nh and vrvr self- 
energy diagrams respectively show a significant effect on the cross-section and consequently 
on the temperature dependence of transport coefficients. The effect of early chemical freeze- 
out in heavy ion collisions is implemented through a temperature dependent pion chemical 
potential which also enhances the medium effect. Results for (/s and rj/s also show a 
significant medium dependence in this scenario. 

The viscous coefficients and their temperature dependence could affect the quantitative 
estimates of signals of heavy ion collisions particularly where hydrodynamic simulations are 
involved. For example, it has been argued in [45] that corrections to the freeze-out distri- 
bution due to bulk viscosity can be significant. As a result the hydrodynamic description 
of the pt spectra and elliptic flow of hadrons could be improved by including a realistic 
temperature dependence of the viscous coefficients. Studies in this direction are in progress. 

VI. APPENDIX A 

The integrals over Bose functions which appear in the definitions of thermodynamic 
quantities like energy density, pressure, entropy etc. of a pion gas can be expressed in terms 
of the functions 5'"(z) where z = m/T. Some of these which appear in various expressions 
in this work are given by 




zT^ 



I 
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Sliz) = dV,p-Ef^\p) (A.l) 



where E = y/jFTm^, f^°\p) = [e^(^'^'^ - 1]-^ and dVp = (Pp/{2-nfE. Using the formula 

oo 

[a — 1]^^ = '^{a^^Y these integrals can be converted to sums over infinite series and can be 

n=l 

oo 

compactly expressed as S'^{z) = e''^^'^k~"Kn{kz), Kn{x) denoting the modified Bessel 

fc=i 

function of order n given by 
or 

, , 2"n!(2n- 1) /"'^ , o 3 , 

Kn{x) = „ ^ / rdrir' - x^y-^e-^ . (A.3) 



{2n)\ X" 

Using now the property 



^ -s:{z) = s:-\z) (A.4) 



9(/i/T) - 

the remaining integrals may be easily obtained. The equilibrium formulae for the number 
density, pressure, energy density and enthalpy density are respectively given by 



n = gj dT,EJ^'\p) = ^z't'SHz) 
P = 9 J dT,p^'\p) = ^z'T'S'.iz 



e = ne 



9 J dT,Elf^{p) = ^z'T^lzSliz) - SUz)] 



H = nh = nzT^fl , (A.5) 

where for a pion gas g = 3. In these appendices we have suppressed the subscript V on m 
and /X for brevity. 



VII. APPENDIX B 

Here we show how the left side of the linearized transport equation (jS]) given by 

p^d,f^'\x,p) = -C[<j>] (B.l) 

can be expressed in terms of thermodynamic forces. We write 9^ = Ufj^u^d^ + A^^d'^ = 
u^D + separating the time derivative D and the gradient V^. Note that D ^ d/dt and 
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V„ — )■ di in the local rest frame. On differentiating /^^^ 



{p-u) 



p ■ u 



, yU , P'^ 



^2 DT + D{^)-^Du, 



+p^ 



p ■ u 

2^2 



T 



/(o)(i + /(o)) ■ 
(B.2) 



The terms DT and -D(^) do not appear in the expression of the thermodynamic forces and 
are to be eliminated using equilibrium laws. From the equation of continuity 



df,m = 



where A^^ = nu^, we get 



Dn = —nd^u^ . 



(B.3) 



Again, contracting the equation for energy-momentum conservation with i.e. 



where T^'^^^^ = n[{e + —)u^u'^ — g^'^—] results in the relation 



De = d.u^" . 

n 



(B.4) 



Further, contracting with the projector A^" , 



yields 



Du^" = — V^P . 
nh 



Making use of the relativistic Gibbs-Duhem relation 



and flRi]) in (IRSD we get 



Using now the expansions 



Dh = TD{j;) + hT-^DT . 



(B.5) 



(B.6) 



Oh = § 

dT 



DT 



de 



DT + 



dh 



T 



9(/i/T) 



T 



(B.7) 
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on the left hand sides of eqs. (1B.4P and (1B.6P results in a coupled set of equations for DT 
and D{fi/T). These are easily solved to arrive at 

dh 



DT 



(P/n) 



T- 



T 



de 




hT-l _ dh 




de 




T- 


dh 






T 




dT 






T. 



fx. 



iP/n) 



irp-l _ dh 



de 




Urp-l dh 




de 




T- 


dh 






T 




dT 






T. 



(B.8) 



We next evaluate the partial derivatives of e and h with respect T and {fi/T) using the 
relations in Appendix-A. We get 



dh 
df 

de 
df 

dh 



.^3 



z 



SI 



q2 qO 

On kDr> 



SI 



cl cO 

s] 

'2 



S' 



CO Cl CO 
^2 _ '^3^3 

S' 



{Sly 



difx/T) 
de 



d{f,/T) 



Tz 
-T 



T 



cO cl cO 
'^3 _ '^3*^2 

C2 CO 
_ ^2^2 



+ Tz 



SI 
SI 



cl cO 
'-'3'-'2 

{s^y 



(B.9) 



Putting these in (IB.SP we get 



TD{^) = [iY~l)h-i"'T]d, 



where 



(B.IO) 

(B.ll) 
(B.12) 

(^27^2)' - {si/siy + 2>z~^s'^si/{siy + 2z-'si/si - z-^ ^^'^^^ 

We now come back to eq. (IB.2p in which we replace DT, -D(^) and Du^ using eqs. (IB.Sp 
and ( IB.SP to get 



7" 



(go/g^^)^ - ('g37'5'2)' + ^z-^sisi/{siy + Z-^g37'g2 
^ " {svsiy - (53751)2 + ?,z-^sisi/{siy + 2^-153751 - z- 

+ (5751)2 - (53751)2 + 3^-150537(51)2 + 2^-153751 - z~'- 

S2IS2 + 5^ "^53/52 — 5053/(52 



7 = 



P 



T 



nh 



T' 



{p.uf{l -i)+p-u {(7" - l)h - i"T} 



1 



3 



+ 



d ■ u = 



/B.14) 
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which on further simphfication and rearrangement gives 
[Qd ■ M + p^A^^(p ■ u - h){T-^d,T - Du,) - {p^,Pu){d^u'')]f^''\l + = -TC\<P] (B.15) 



where the abbreviation 



Q = -Im' + (p ■ ufi^ - 7'} + P ■ «{(7" - l)h - 7'"T} 



(B.16) 



and 



3 

p2 



^ 12 

{p^lPu) = Pf,Pu - y^A*!^ - P ■ '^iPtiUi, + Pi,u^ - -gf,uP ■u) + ^{p- ufu^u^ . (B.17) 



VIII. APPENDIX C 



Here we briefly describe how the leading order expressions for the viscosities are obtained 
from the i nteg ral equations f|T3|) and (1141) . Let us first consider the bulk viscosity Q. Fol- 
lowing 2l|, |28| we multiply both sides of eq. ( !T3l) by Laguerre polynomial of order 1/2 and 
degree n and integrate over dVp to get 



(C.l) 



where r = (3{p ■ u — m) and 



\;ldT,f^'\p){l + f^'\p)}QL]!^i 



riT 



(C.2) 



and the abreviation 



[F,G] = ^J dT,dT,dT,, dT,,f^^\p)f^^\k){l + f^^\p')}{l + f'\k')}6{F)6{G) W (C.3) 



with 



5(F) = F{p) + F{k) + F{p') + F{k') 
5(G) = Gip) + Gik) + Gip') + Gik') 



(C.4) 



Expanding A as 



00 



m=0 



(C.5) 
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and putting in ( IC.ip the expression for bulk viscosity to first order is obtained as 



C = T 



^22 



(C.6) 



where 02 is given by eq. (1C.2|) with n = 2. The denominator 022 in ( 12T|) can be expressed 
using (IC.SP as 

a22 = [Lf (r),Lf (r)] (C.7) 

and involves a 12- dimensional integral. This is considerably simplified by making a change 
of variables as described in 2l|, |28|] reducing the number of integrations essentially to five so 



that 



m 



0-22 



dx / dtp 



de 



2tt 



dQ 



1 da 



X 6{Lr{r)}6{Lr{r)}-^{^,e) 



X sinh^ X sinh^ cosh'^ sin G sin 9 

Using now 

f^'\p)f^'\k){l + f^'\p')}{l + f'\k')} 



(C.8) 



2fi/T p2z cosh -)/) cosh x 



(eS-l)(e^-l)(eG-l)(e^-l) 
where the exponents in the Bose functions are given by 

E = z(cosh -i/; cosh X — sinh-^sinhx cos 6') — /i/T 

F = z(cosh t/' cosh X — sinh'i/'sinhx cos 6'') — /i/T 

G = E + 2z sinh ip sinh x cos 6 

H = F + 2z sinh sinh x cos 6' 



and 



we finally get 



6{lI^\t)}6{L2'\t)} = ;2^(sinh^sinhx)^(cos^^ - cos'^') 



1/2 



2 /i/\2 



022 = 22:^/3 (z) 



(C.9) 



(C.IO) 



(C.ll) 



The relative angle 9' is defined by, cos 6*' = cos 6' cos 6 — sin 6* sin 9 cos . Also, on simplifi- 
cation, the integral 02 is finally given by 



«2 = ---^) + (| + 7|){(i-7")| + 7"V^} 

o O2 O2 



(C.12) 
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The calculation for the shear viscosity r] follows a similar procedure and is more involved 
because of the tensorial nature of C^j, = C{p^piy). In this case eq. f lT^ is multiplied by 
Laguerre polynomial of order 5/2 on both sides and then integrated over dTp to get the 
formal expression 

[C{r){p,p.), (r)(p>0] = -Tin (C.13) 



m 
n 

where 



In 



■ / rfr,/(°)(p){l + /(°)(p)}L^/^(r)(p,p.)(p'^p^) . (C.14) 



mnT^ 
Writing in this case 

oo 

C{r) = Y: CmL'l\r) (C.15) 

m=0 

and putting in fIC.lSp . the first approximation to rj is obtained as 



10 Coo 

where 



T 

^=T7T— (C.16) 



1 



coo = ^[^f (r)(p,P.),Lf (r)(pV)] . (C.17) 



Simplification along the lines of [2l|, |28| finally yields 



16 

Coo = 16/i(^) + IQhiz) + y/3(^) . (C.18) 



The integrals Ia{z) are defined as 

/ dijj cosh^ ■^Z' sinh"^ T/^ / (i6 sin 6- — (■?/', 9) / dcj) 
Jo Jo 2 ail Jo 

POD riT 2z cosh V cosh X 



/ f^) = f e(-2M/T) 



where the functions Mn stand for 



Ml (^,6) = l-cos^0 , 
Ms (^,6) = cos2^ + cos2^'-2cos^cos0'cose , 
M3(^,0) = [cos^ - cos^ (C.20) 
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